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The interplay between quantum-mechanical properties, such as coherence, and classical notions, 
such as energy, is a subtle topic at the forefront of quantum thermodynamics. The traditional Carnot 
argument limits the conversion of heat to work; here we critically assess the problem of converting 
coherence to work. Through a careful account of all resources involved in the thermodynamic 
transformations within a fully quantum-mechanical treatment, we show that there exist thermal 
machines extracting work from coherence arbitrarily well. Such machines only need to act on 
individual copies of a state and can be reused. On the other hand, we show that for any thermal 
machine with finite resources not all the coherence of a state can be extracted as work. However, 
even bounded thermal machines can be reused infinitely many times in the process of work extraction 
from coherence. 
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I. INTRODUCTION 

Scientia potentia est, knowledge is power, the Latin 
aphorism goes. This could not be more true in thermo¬ 
dynamics, where knowledge about the state of a system 
can be exploited to our advantage to extract work from it 
[1, 2]. In quantum mechanics states of maximal knowl¬ 
edge are called pure states. A peculiar feature of the 
quantum world is that, due to the superposition princi¬ 
ple, even for such states there are many questions that 
cannot be answered sharply. In thermodynamics we are 
especially interested in energetic considerations and so 
an odd place is taken by pure states that are a super¬ 
position of different energy states. This is because, de¬ 
spite the fact that we possess full knowledge about the 
system, our possibility of predicting the outcome of an 
energy measurement can be very limited. 

In standard quantum-mechanical considerations this is 
not a issue, because we can always reversibly transform 
a pure state into any other pure state by unitary dy¬ 
namics. A basic task of thermodynamics, though, is the 
book-keeping of all energy flows from and out of the sys¬ 
tem, and there is no reversible transformation mapping 
a superposition of different energy states into an eigen¬ 
state while strictly conserving energy. Hence, we are left 
to wonder whether the “scientia” of having a pure state 
with quantum coherence can be converted into “poten¬ 
tia” of extracted work, while being limited by the law of 
energy conservation. 

More precisely, we analyze work extraction from quan¬ 
tum coherence^, in the context of the theory of thermody¬ 
namics of individual quantum systems, currently under 
development [3-16]. The aim of the theory is to provide 
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^ Here, and in the rest of the paper, we use the term “coherence” in 
the sense of “superposition of states belonging to different energy 
eigenspaces”. 


a suitable theoretical framework for our increasing abil- 
ity to manipulate micro- and nanoscale systems [17-21]. 
Such general framework could also help approaching re¬ 
lated questions, such as the role of quantum effects in bi¬ 
ological systems [22-24] and the link between thermody¬ 
namics and quantum information processing. Although 
this field has recently seen a great number of contribu¬ 
tions, most of the previously mentioned works do not 
incorporate the possibility of processing a state in a su¬ 
perposition of energy eigenstates. It was only relatively 
recently that the role of coherence in thermodynamics 
has been looked at more closely [25-34]. 

In this paper we first set the scene by presenting 
existing approaches to the problem of work extraction 
from the coherence of quantum systems. We argue that, 
within the regime of individually processed systems, the 
current approaches fail to account for all the resources 
used during the work extraction protocol. The typical 
assumption is to use a classical external field that expe¬ 
riences no back-reaction [35, 36]. However, this does not 
allow for a full accounting of the thermodynamic cost of 
maintaining the field. Although this cost may be small 
in a single use, it has to be accounted for since the work 
gain will also be small. Hence we propose an alterna¬ 
tive framework that aims for a careful book-keeping of 
resources. 

In particular, we use the notion of a thermal machine 
[37] , a device of bounded resources that can be used to 
manipulate thermodynamical systems and perform tasks 
such as work extraction. Our thermal machine incorpo¬ 
rates the use of an ancillary system carrying coherence 
(henceforth called reference system), introduced into the 
context of coherence manipulation in thermodynamics in 
[26]. It also includes a battery system where work can 
be stored, or transferred to the reference when necessary. 
A crucial question addressed in this work will be how to 
use the thermal machine in a repeatable way, i.e. with¬ 
out deteriorating it. We make use of an important result 
of Johan Aberg, showing that reference systems can be 
used repeatedly to manipulate coherence [31]. However 
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the reference needs to be repumped for the machine to 
continue to operate. The work cost of repumping needs 
to be taken into account, which can mean that the ther¬ 
mal machine will not be able to extract all the available 
work from coherence. Nonetheless, we will find that we 
can come arbitrarily close, by choosing the amount of co¬ 
herence resources carried by the reference system in the 
machine appropriately large. 

On the other hand, for any given thermal machine, 
we will prove that one can never extract all the avail¬ 
able work. We will show that coherences of individual 
quantum systems can be exploited to enhance the per¬ 
formance of work extraction protocols (both in the aver¬ 
age and the single-shot sense), but not to the extent that 
could be expected in the “classical” limit. Moreover, the 
work extraction protocol we provide does not deteriorate 
the thermal machine. 


II. COHERENCE AND WORK 
A. Setting the scene 

Let us start by introducing the framework that we will 
use throughout this paper and collect the core assump¬ 
tions that our results rest upon. We want to study the al¬ 
lowed thermodynamic transformations by explicitly mod¬ 
elling any coherence resources being used. There are two 
ways in which coherence can enter the thermodynamics 
of the systems under consideration. This can happen ei¬ 
ther explicitly, by transferring it from an external system 
with quantum coherence (a trivial example being a swap 
operation between the system and an ancillary coherent 
state); or implicitly, by allowing operations that do not 
conserve energy (e.g., |0) —> |-|-), where the Hamiltonian 
is given by Hg = |1)(1|) or conserve it only on average 
(e.g., |1) ^ (|0) + |2))/y2, with Hs = |1)(1| + 2 |2)(2|). 
Therefore, we will only allow for those transformations 
that do not implicitly introduce coherence: 

Assumption 1 (Allowed transformations). The set of 
allowed transformations is given by all (strictly) energy¬ 
preserving unitaries, i.e., unitaries that commute with 
the total free Hamiltonian of the system. The use of all 
ancillary systems should be explicitly accounted for.^ 

We will also take a closer look at an alternative approach 
in Sec. HD 1 and explain why we find it not satisfactory 
for the aims of the present work. 

In this paper we focus on the task of work extraction 
from quantum systems with coherence. We do not aim 
here to settle the long-standing issue of what is an ap¬ 
propriate definition of work in quantum thermodynamics 


^ As we will see, thermal ancillas are the only ones that can be 
freely introduced without trivialising the problem of work ex¬ 
traction. 


(see, e.g., [38, 39]). For the scope of this paper we will as¬ 
sume, for the sake of simplicity, that the following holds 
for classical (incoherent) states: 

Assumption 2 (Average work, incoherent states). Let 
PS be a quantum state of the system described by Hamil¬ 
tonian Hg, with pg being incoherent in the energy eigen- 
basis. Then, in the presence of a heat bath at tempera¬ 
ture T, an average amount of work {W){ps) equal to the 
change of free energy of a state can be extracted from it: 

{W){pg) = AF{pg):=F{pg)-F{jg), (1) 

where F{a) := Tr {crFfg) — kTS{a), S{-) is the von Neu¬ 
mann entropy and jg = jZg is a thermal state 

with Zg being the partition function of the system. 

This formula, consistent with traditional thermodynam¬ 
ics, has been obtained using work extraction models that 
differ in details, but agree on the result [9, 35]. For ex¬ 
ample, in Ref [9] the work extraction protocol is based 
on two elementary processes: level transformations (that 
change the eigenvalues of the system Hamiltonian Hg) 
and full thermalisation with respect to the current sys¬ 
tem Hamiltonian (through thermal contact with a bath 
at temperature T). Average work is then defined as the 
average change in energy during level transformations 
(the “unitary” steps) and, if initial and final Hamilto¬ 
nian coincide, Eq. (1) is recovered. Here we focus on 
the problem of extending Eq. (1) to quantum states with 
coherence. The results of the present paper apply to any 
definition of work satisfying Assumption 2. 

The problem of work extraction can also be studied 
in the so-called single-shot regime. This means that one 
is interested in single instances of the work extraction 
protocol, instead of average quantities. To explain this 
more precisely, let us refer again to the model introduced 
in Ref. [9] that we have summarised above. Extracted 
work can then be seen as a random variable, maximis¬ 
ing the average of which yields Eq. (1). However, we 
may instead ask what is the maximum amount of de¬ 
terministic (i.e., fluctuation-free) work that can be ex¬ 
tracted during a single instance of the protocol, while 
allowing the failure probability e. In Ref [9] it was 
shown that for incoherent states this quantity is given by 
FoiPs) - Fijg), where F^{pg) = -kTlog Z^ is a single¬ 
shot free energy defined as follows. Given a subset A of 
the indices {i} labelling the energy levels of the system, 
define Z{A) = where Ei are the eigenvalues 

of Hg. Then = minA{Z(A) : ft > 1 - e}. where 

Pi = {Ei\pg\Ei). This result is also in agreement with 
other work extraction models based on thermal opera¬ 
tions [7]. Hence, for the single-shot scenario we can use 
the following assumption: 

Assumption 3 (Single-shot work, incoherent states). 
Let pg be a quantum state of the system described by 
Hamiltonian Hg, with pg being incoherent in the energy 
eigenbasis. Then in the presence of a heat bath at tem¬ 
perature T, using a single-shot protocol one can extract 
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a sharp amount of work with failure probability e: 

^tsiPs) = AFJ(ps) := Fq{ps) - F{js)- (2) 

Once again, our aim is to extend Eq. (2) to quantum 
states with coherence and our results apply to any defi¬ 
nition of single-shot work satisfying Assumption 3. For 
the sake of brevity in the remaining of this paper we will 
only write “extracting work equal to the free energy”, 
omitting “in the presence of a heat bath at temperature 
T”; however, this is how our claims should be understood. 

In thermodynamic considerations thermal Gibbs states 
are the only ancillary states that can be introduced with¬ 
out the need for careful accounting. In fact, one can 
show that using energy-preserving unitaries (in accor¬ 
dance with Assumption 1), a thermal state is the only one 
that can be introduced for free without allowing the pro¬ 
duction of every incoherent state [8]. Clearly, if this was 
possible, then from Assumptions 2 and 3 one could ex¬ 
tract infinite amount of work, thus trivialising the theory. 
Hence, the most general thermodynamic transformations 
that can be performed without using extra resources are 
given by: 

1. adding a bath system in a thermal state 'Je with 
arbitrary Hamiltonian He and fixed inverse tem¬ 
perature P = l/fcT, 

Ps^ PS®1E, 7£; = e"^^®/Tr ; (3) 

2. performing any global unitary that conserves to¬ 
tal energy, i.e., that commutes with the total free 
Hamiltonian of the system and baths, in accordance 
with Assumption 1; 

3. discarding any subsystem. 

The set of quantum maps acting on a system that arise 
from combining the transformations described above is 
known under the name of thermal operations [3, 26]. 

B. Work-locking 

The aim of this work is to begin with a system initially 
in a state with coherence ps, and finish with a thermal 
state 7 s, while optimally increasing the free energy of a 
battery (storage) system. The initial and final battery 
states, pb and should be incoherent, so that using 
Assumptions 2 and 3 we can achieve the coherence to 
work conversion that we are looking for. Schematically: 

PS®Pb^1s®Pb- (4) 

Without the use of an ancillary resource state the above 
transformation is given by a thermal operation. Note 
that thermal operations commute with the dephasing 
channel T) [26] that removes all coherence from a quan¬ 
tum state, 

V[<7) :=^Tr(n,a)n„ 


where H^ are the projectors on the energy eigenspaces of 
the system under consideration. Hence, we get that if 
the transformation described by Eq. (4) is possible, also 
the following one is: 

V{ps)®PB -^1 s®p'b- 

This implies AF{pb) < AF{'D{ps)), because F is non¬ 
increasing under thermal operations. From Assump¬ 
tion 2 we then have {W){ps) < {W){V{ps)). A similar 
argument gives also W^^{ps) < Wgs{'D{ps)); note that 
in both cases the bound is achievable because dephas¬ 
ing is a thermal operation. This phenomenon was ob¬ 
served before [7, 10] and was called “work-locking” in 
[28]. Work-locking highlights that, despite contributing 
to the free energy of the state, quantum coherence does 
not contribute to work extraction: it is “locked”. It also 
shows, in agreement with [31], that the standard formula 
{W){ps) = AF{ps) applied to every state (also the ones 
with coherence), implicitly assumes the access to an ex¬ 
ternal source of coherence. In this paper we revise the 
problem of extracting work from coherence, clarifying the 
role of this external source of coherence. To summarise 

Central Question. To what extent can Eq. (1) and 
Eq. (2) be extended to arbitrary quantum states with 
coherence, while explicitly accounting for coherence re¬ 
sources (ancillary systems)? 

C. Different thermodynamic regimes 

With an increasing interest in the thermodynamics of 
non-equilibrium quantum systems, an important distinc¬ 
tion to make is between “single-shot” statements, which 
are valid for every run of the protocol, and “many-runs” 
statements, valid in the case of a large number M of 
runs. In the asymptotic regime M —> oo one is focused 
on studying average quantities (like average extracted 
work), which is justified by the fact that the fluctuations 
around the average can be made negligible in the limit 
of a large number of runs of the protocol (which is often 
the situation of interest in the study of heat engines). 
On the other hand, although the expected amount of ex¬ 
tracted work can be studied in a single-shot regime [35], 
it potentially carries little information about the system 
at hand due to the large fluctuations of non-equilibrium 
thermodynamics [9]. Instead, the focus in the single¬ 
shot regime is typically on probabilistic work extraction 
protocols that guarantee precise and sharp amount of 
work with a finite probability of success or some mini¬ 
mum amount of guaranteed work [5, 7, 9, 11]. On top of 
this classification we can also differentiate between “indi¬ 
vidual processing” scenarios, in which a single (possibly 
nanoscale) system undergoes a thermodynamic process 
on its own; and “collective” scenarios, in which TV > 1 
copies of a state are processed together (the N ^ oo 
limit is considered in [26]). 

This classification of the thermodynamic regimes in 
which work extraction can be analysed is presented in 
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FIG. 1. Thermodynamic regimes. Work extraction 
protocols can be investigated in different thermodynamic 
regimes. These can be classified by the number of systems 
that are processed at each run of the protocol (individual vs 
collective) and the number of times the protocol is repeated 
(single-shot vs many-runs). The green background indicates 
that in a given regime the maximal amount of work that can 
be extracted is consistent with traditional thermodynamics. 

Fig. 1. The work-locking described in the previous sec¬ 
tion is a feature appearing in the regime of individually 
processed quantum systems. On the other hand, allowing 
for collective processing of the systems or for an ancillary 
quantum memory, one can “unlock” work from quan¬ 
tum coherence (see Appendix A for a short description 
of this point). For arbitrary collective processes, one only 
needs to use a sublinear amount of coherence in the refer¬ 
ence system, meaning that its consumption does not con¬ 
tribute to the average work consumed or produced in the 
N ^ oo limit [26]. In this limit coherence plays no role 
as, e.g., for N identical qubits F{p^^) « F{'D{p®^)), 
with a deficit scaling as \ogN/N [28]. Hence, in this pa¬ 
per we are interested in the thermodynamics of individual 
quantum systems {N = 1, the upper half of Fig. 1). 

D. Individual processing regime 

1. Average energy conservation 

In [35] sharp energy conservation, as expressed by the 
unitary dynamics commuting with the total Hamiltonian 
FI tot, was replaced with the condition that such dynam¬ 
ics only keeps the first moment (iFtot) constant. Under 
this weaker condition it was shown that an amount of 
energy equal to the free energy difference AF{ps) can 
be extracted on average from a system in an arbitrary 
quantum state ps- 

The elegance and appeal of this is that it recovers a 
clear thermodynamic meaning for the free energy of an 
individual quantum system. However, several problem¬ 
atic issues can be raised. Firstly, if one is interested in 
analysing the class of allowed quantum operations, then 
in the average-energy scenario this set depends on the 


particular state one is processing, which is conceptu¬ 
ally less appealing and technically problematic from a 
resource-theoretic perspective. 

Secondly, restricting energy considerations to the 
first moment analysis can hide arbitrarily large en¬ 
ergy fluctuations described by higher moments, that are 
not explicitly modelled, but may be highly relevant. 
To see this consider a unitary Have mapping a state 
IV^ 02 ) := (jo) + \2))/t/2 to Jl), which preserves energy on 
average (here jn) is the energy eigenstate corresponding 
to energy n). Since microscopically all processes are ulti¬ 
mately energy-conserving, C/ave must be realized through 
a joint energy-preserving unitary U involving IV’ 02 ) and 
some ancillary state pA, e.g., the state of the battery, 

U{\ll)02){'>Pt)2\® ® p'a- 

In any such process the energy fluctuations of the ancil¬ 
lary system must increase. Specifically, denoting by F[ 
the Shannon entropy of the outcomes of an energy mea¬ 
surement one gets a strict inequality: FI{p'j^) > FI {pa) 
(see Appendix B for details). As an example consider the 
ancillary system prepared in the energy eigenstate jm), 
so that an energy measurement would give a sharp out¬ 
come. Then, while the system is transformed from 1'0O2) 
into jl), the ancilla must be transformed into a superposi¬ 
tion of energy eigenstates (Jm -I- 1) -I- [m — l))/-\/2. Hence 
an energy measurement would show fluctuations in the 
final state of the ancilla. It is important to note that the 
protocol that extracts work from coherence within the 
framework of average energy conservation [35], necessar¬ 
ily creates such extra fluctuations, however these are not 
explicitly modelled within the formalism used. As we will 
see it is exactly due to these fluctuations that our proto¬ 
cols require work to be invested in restoring the ancillary 
state. 

Finally, as the fluctuations created by operations that 
conserve energy only on average remain outside the for¬ 
malism, one cannot properly account for the fluctuations 
in the extracted work outside the asymptotic regime. 

2. Repeatable use of coherence resources 

As already mentioned in Sec. HB, in the presence 
of energy conservation and without additional coherence 
resources, work-locking prevents us from extracting work 
from the coherence of individual quantum systems. One 
could then stay within the framework of strict energy 
conservation, but allow for the use of an extra source of 
coherence. We refer to this extra system as the reference. 

At one extreme one could allow for the use of an infinite 
source of coherence (an “unbounded” reference frame 
[40, 41])^, that entirely negates the constraints and ex¬ 
periences no back-reaction from its use on the quantum 


® Not to be confused with a reference described by a Hamiltonian 
unbounded from below, which is unphysical. 
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system. As suggested in [28], in such case we should 
be able to extract all the work from coherence. How¬ 
ever, one might worry that this involves the accounting 
“oo — c = oo” , with c being some finite resource con¬ 
sumed from an infinitely large reference system. Indeed, 
the use of such an unbounded reference allows us to simu¬ 
late the operations from the previous section (conserving 
energy only on average) [41], and hide the arising extra 
fluctuations in the infinitely big reference system. This 
semiclassical treatment is typical for many standard ap¬ 
proaches that assume the existence of a classical field 
experiencing no back-reaction from the system [35, 36] 
and works well in many circumstances. However we are 
interested in the regime in which the thermal machine 
itself may be a microscopic quantum system. Hence, it 
seems more reasonable to firstly consider the reference 
as a quantum system with finite coherence resources - 
a “bounded” reference frame - and only then study the 
limit of an unbounded reference (recent works in this 
spirit and the discussion of semiclassical approaches can 
be found, e.g., in [29, 38, 42]). 

Definition 1 (Reference). We consider a refer¬ 
ence (or coherence reservoir) given by an infinite¬ 
dimensional ladder system described by Hamiltonian 
Hr = J2^=o We characterise the state pn of the 

reference through two numbers, ((A),M). The hrst pa¬ 
rameter, (A), measures the coherence properties of the 
reference and is given by 

(A)=Tr(pflA), A = (A + At)/2, (5) 

where A is the shift operator A = 

We have that (A) < 1 and the limit case (A) = 1 is 
called unbounded or classical reference. The second pa¬ 
rameter, M, describes the lowest occupied energy state, 
M = minjn : {n\pfi\n) > 0}. 

Examples of a sequence of references that come arbitrar¬ 
ily close to a classical one are uniform superpositions of 
L energy states when L —> oo or coherent states with 
arbitrarily large amplitude. The use of (A) and M as 
relevant quality parameters will soon become clear. 

Results from the field of quantum reference frames 
[40, 41, 43-45] suggest that the back-reaction experienced 
by the reference will necessarily deteriorate it and con¬ 
sume the resources. However, if the usefulness of the 
reference or field is continually degraded during the work 
extraction process, we cannot claim that we are present¬ 
ing a protocol performing work extraction from the state 
alone, as extra resources are consumed. Similar prob¬ 
lems arise if free energy is continually taken away from 
the reference. 

In this paper we propose the following approach. We 
allow for the use of additional coherence resources as part 
of our thermal machine, but demand that they are used 
repeatably in the following sense: the performance of our 
reference-assisted protocol, while operating individually 
on the n-th copy of the system, must be the same as while 


operating on the {n + l)-th copy, for all n S N. In other 
words, repeatability means that the reference’s ability 
to perform the protocol never degrades, but crucially its 
state is allowed to change. Essentially this means that 
despite that the free energy of the reference can fluctu¬ 
ate and its coherence properties change, it can be used 
indefinitely to repeat the same protocol. To design such 
a protocol we employ the recent surprising result of [31] 
that shows how a coherence resource can be used repeat¬ 
ably to lift the symmetry constraints imposed by energy 
conservation.^ However, as we will see the protocol in 
[31] requires continuous injection of energy into the ref¬ 
erence (we do not allow the Hamiltonian of the reference 
to be unbounded from below, as in [47]). Hence, it is not 
immediately obvious that net thermodynamic work can 
be extracted from coherence. 

In what follows we introduce a general protocol that 
processes quantum systems individually and allow us to 
extract work from their coherence. We then focus on two 
variations of it. The first one can come arbitrarily close 
to extracting all the coherence as average work with ar¬ 
bitrarily small failure probability, provided we make the 
coherence resources of the reference system in the ther¬ 
mal machine large enough. However, if one does not 
have access to arbitrarily large coherence resources, this 
variation of the protocol does not guarantee perfect re¬ 
peatability. Therefore, we examine a second variation 
that is perfectly repeatable even for bounded references. 
We show then that the performance of work extraction in 
both the single-shot and asymptotic regimes is enhanced 
only if the quality of the reference (defined further in the 
text) is above a certain threshold. 

III. THE PROTOCOL 

We analyze work extraction from pure qubit states 
with coherence, 

lV>) = ^r^|0) + v^e-*^lI), pe(0,I). (6) 

Without loss of generality we can set Hs = |1)(1| and 
ip = 0 (rotations about the z axis of the Bloch sphere 
conserve energy). Our aim is to unlock work from co¬ 
herence through the repeatable use of a thermal machine 
containing a reference, while processing each copy of lip) 
individually. In Table I the results we obtain within this 
framework are schematically compared with the ones ob¬ 
tained within the frameworks presented in the previous 


The work [31] actually uses the word “catalysis”, but we prefer 
to use the word repeatability/repeatable to avoid suggesting that 
there is no change in the state of the reference. Recall that 
traditionally a catalyst is a system in a state y that enables 
p 0 X 0 X 5 despite p ^ <7 being impossible (see, e.g., [ 8 , 

46]). Repeatability, on the other hand, only requires the auxiliary 
system to be as useful at the end as it was at the beginning, while 
its state may change. 
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Single-shot 

Asymptotic 

Average energy 

conservation 

{W) = AF [35] 
Large fluctuations [9] 

{W) = AF [35] 

Strict energy 

conservation 

e(l7)) = e( 7 s), {W) = 0 [7, 10, 28] 

Strict energy 
conservation w/ 
resource used 

repeatably 

e(| 7 )) = 0 

unbounded reference 

e(l7)) < c(7s) 
bounded reference 

(IT) = AF 
unbounded reference 

(IT) < AF 
bounded reference 


TABLE I. Individual processing protocols extracting 
work from I 7 ). (W) denotes the average work that can be 
extracted from the coherent thermal state I 7 ) and e denotes 
the error probability of a single-shot work extraction from 
a given state. A thermal state of the system is denoted by 
7 s. Note that under operations strictly conserving energy, 
no work extraction protocol on I 7 ) can outperform a work 
extraction protocol on 7 s, as the two states are indistinguish¬ 
able. 


section for the paradigmatic example of a qubit in a “co¬ 
herent Gibbs state” I 7 ) given by: 

I 7 ) = \/l - rlO) -b V^|l), (7) 

with (1 — r, r) being the thermal distribution for the sys¬ 
tem, so that 27(1 7 X 7 1) = 7s. 

The extraction of non-zero work from \'il)) then requires 
a thermal machine containing a reference state pn and 
implementing an energy-conserving unitary V : 

p'sR = PR)y\ 

satisfying the following: 

1. The system is pre-processed to a new state 
p'g = Trjj {p'sr) that allows for better work extrac¬ 
tion than from the initial state |' 0 )- 

2. The final reference state p’j^ = Trs {p’gj^ can be 
processed into a state p'^ (perhaps using some of 
the extracted work) in such a way that the repeata¬ 
bility requirement is satisfied. 

3. No collective operations, at any stage of the pro¬ 
tocol, are allowed on multiple copies of |X’) and no 
quantum memory (in the sense of Appendix A) is 
used. 



FIG. 2. The basic protocol. The evolution of the system 
from the initial state \'ip) to the final state is depicted on the 
Bloch ball, in blue and red respectively. The evolution of the 
reference from the initial state (smaller blue blob) throughout 
the protocol (red blobs) is depicted on the energy level ladder. 


the reference pR. The unitary acting on the joint 
system SR is chosen as in [31] to be: 


y(;7) = |o)(o|G|o)(o| + ^yXi7), (s) 


with 


1 

V,{U)= Y. (n| U \m) \n){m\ 0 X — n){l — m\. 

n,m=0 


We choose 


U = 


y/p -V^-p\ 

y/^-p y/P ) ’ 


so that U rotates the qubit system from \'ijj) to |1).® 

2. Work extraction. The system is now in a state 
p’g and, due to work-locking (see Sec. II B), is indis¬ 
tinguishable from its dephased version in any work 
extraction protocol. So without loss of generality 
we can use the dephased version 


V{p's) = {l-q)\Om+q\l){l\. (9) 


Now, depending on the considered regime, single¬ 
shot or average (asymptotic) work can be extracted 
from the dephased state and stored in the thermal 
machine. 


3. Repumping. The back-reaction changes the state 
of the reference into p'^. Using part of the extracted 


A. The explicit work-extraction protocol 

A protocol satisfying the introduced requirements con¬ 
sists of the following steps (see Fig. 2): 

1. Pre-processing. The system jX:) interacts 
through an energy-preserving unitary V(U) with 


® This interaction corresponds to a modified Jaynes-Cummings 
model (with excitation-dependent coupling strengths). However, 
it can also be approximately realized within the standard Jaynes- 
Cummings model using a reference in a coherent state 10}, with 
|a| large enough (for details see Supplementary Material Sec. V 
in [31]). 
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work (stored in the battery during the previous 
step) we can repump the reference to shift it up: 

p'r p'k '■= ^p'n^k (10) 

with the details of how to perform such operation 
given in Appendix C. We will describe how often 
we perform the repumping while analysing different 
variations of the protocol. 

4. We can repeat the protocol using p'^ and a fresh 
copy of IV^). 

B. Performance 

During the pre-processing stage the joint unitary V(U) 
approximately induces U on the system: 

pk^u{mi;\)u^ = m\. 

The degree to which the above equation holds depends 
on the quality of the reference as defined in [31]. In par¬ 
ticular, the system final occupation in the excited state 
9 = (l|Ps |1) is given by 

9 = 1- 2p(l - p)(l - (A)) - (1 - p)^Roo, (11) 

where i?oo = (0|Pfl|0) and A is the quality parameter 
from Definition 1. From Eq. (11) it is easy to see that 
9 —>■ 1 when Rqq 0 (i.e., when the reference quality 
parameter M > 0) and (A) —>• 1. Therefore, Rqq and 
(A) are operationally well-defined quality parameters of 
the reference, because they directly measure the ability 
of the reference to induce the unitary U that we want 
to perform on the qubit. At the same time the reference 
undergoes a back-reaction induced by the joint unitary 
V{U). This is described by the following Kraus opera¬ 
tors: 

^0 = (1 - -P|0)(0| + \/p(l -p)(II - A), (12a) 

Ai = pi -I- (1 - p)A'^, (12b) 

so that the reference state after performing the pre¬ 
processing stage is given by: 

p'r = ^oP^Aq -I- AipuA^. (13) 

From Eqs. (9) and (11) the only two parameters rel¬ 
evant for work extraction are the reference population 
in the ground-state, i?oo, and the parameter (A). Using 
the Kraus operators specified by Eqs. (12a)-(12b), the 
change in (A) during the pre-processing stage (i.e., the 
difference between the final and initial value of (A)) can 
be computed and is found to be: 

(5(A) = (1 - p) [(1 -p - 2p)Re(Aoi) - y/pRm ], (14) 

where as before = {i\pR\j)- A sufficient condition for 
(A) to stay constant is i?oo = 0; i-O., M > 0. Therefore, 
if the initial state satisfies Rqo = 0 exactly, performing 


the pre-processing stage does not change (A), as noted 
in [31]. We require step (3) of the protocol to ensure 
that (Ojp^jO) = 0. If this is the case, at the end of the 
protocol we are left with the reference described by the 
same quality parameters Rqo and (A) as at the beginning, 
and the reference p'^ is as good as pR within the protocol. 

Finally, because the state of the reference changes, its 
free energy can fluctuate. However, notice firstly that 
the reference has Hamiltonian bounded from below, so 
for fixed average energy it has a finite amount of free en¬ 
ergy. Secondly, repeatability requires that the reference 
can be used an arbitrary number of times and the per¬ 
formance of the protocol never changes. It is then easy 
to see that on average the free energy of the reference 
cannot be extracted as work, as this would be incom¬ 
patible with repeatability. It can be shown that in the 
worst-case scenario the free energy change in the refer¬ 
ence fluctuates around zero (see Appendix D for more 
details). Moreover, these fluctuations vanish in the limit 
of an unbounded reference (A) —>■ 1, as then the entropy 
of the reference stays constant, while its average energy 
increases. Therefore its free energy must increase at ev¬ 
ery step of the protocol. 

IV. FUNDAMENTAL LIMITATIONS OF 
COHERENCE TO WORK CONVERSION 

How well does the above approach do in terms of work 
extraction? Here, we first emphasize the limitations on 
work extraction from coherence that arise due to the ref¬ 
erence being bounded, i.e., when we have access to lim¬ 
ited coherence resources. More precisely, we will explain 
why the use of a bounded reference does not allow one to 
extract from a state with coherence the average amount 
of work equal to the free energy difference. 

However, we will then also show that one can construct 
a series of bounded reference states that come arbitrarily 
close to extract the free energy AF(|^)), with protocols 
arbitrarily close to perfect repeatability®. Thus, we will 
prove that in the limit of an unbounded reference all co¬ 
herence can be converted into work in a repeatable way. 
The limit case does not generate any entropy in the ref¬ 
erence system and, being a reversible transformation, is 
optimal. 

A. Limitations of bounded thermal machines 

In order to illustrate the limitations arising from us¬ 
ing a bounded reference we will consider a particular 
model of work extraction from coherence described in 


® A similar result appears in [31], however it was based on using a 
reference system described by a doubly-infinite ladder Hamilto¬ 
nian. This left open the question if this limit is achievable by a 
system with a physically realisable Hamiltonian. 





Refs. [35, 36]. It has been proved there that the free 
energy difference AF{ps) can be extracted from a sys¬ 
tem ps as work if one allows for the use of operations 
that conserve the energy only on average. Let us briefly 
recall the protocol used to achieve this. It is com¬ 
posed of two stages; first, given a state ps, work is ex¬ 
tracted from coherence. The resultant state is given by 
T^{ps) (recall that V denotes the dephasing operation 
in the energy eigenbasis). Second, work is extracted 
from the incoherent state V^ps). In accordance with 
Assumption 2, the latter process extracts AF{'D{ps)). 
Hence, the extraction of the full free energy AF{ps) 
from a state ps is equivalent to the possibility of ex¬ 
tracting F{ps) — F{V{ps)) from the coherence of ps- 
Notice that this quantity coincides with kTA{ps), where 
^{ps) = S{ps\\F>{ps)) is a known measure of quantum 
coherence [48], and it quantifies the amount of free en¬ 
ergy stored in coherence [28]. Hence, the amount of work 
that needs to be extracted on average from the coherence 
of a quantum state to achieve the free energy extraction 
limit for arbitrary quantum states is 

W'coh(ps) = kTAips). (15) 

Without loss of generality we can write any state ps as 
J2nPn IV^n)(V’ra| with pn +1 < Pn- Let US also denote the 
Hamiltonian of the system by FIs = \En){En\- In 

the protocol that allows to extract work from coherence 
in Ref. [35] , the system ps interacts with a weight system 
in a gravitational field via the unitary 

Ua-ve — E \En){^n\<S>Ts^, (16) 

n 

where Lg^ is the shift operator on the weight system that 
shifts it in energy by e„ = {ipnWsHn) - E^- 

As Uave does not strictly conserve energy, by Assump¬ 
tion 1 it is not a free thermodynamic operation. One can 
instead ask if it can be achieved by an energy-preserving 
unitary R(C/ave) on a larger system that exploits some 
ancillary system pA (that includes a battery). In other 
words, we are looking for the energy-preserving unitary 
y (t^ave) such that 

£{ps) ■= Tr^l {V{Ua.ve){pS ® PAWiUaveY) = Ui^vePS U^^e- 

(17) 

It is easy to show that due to the imposed constraints, 
the ancillary system pA must carry quantum coherence. 
In fact, if Pa were incoherent, then the left-hand side of 
Eq. (17) would be a time-translation covariant quantum 
map (meaning that [f,I?] = 0) [33], whereas the right- 
hand side is not. 

Now the crucial point is that Eq. (17) cannot hold 
exactly unless pA contains an unbounded reference. If 
PA is bounded, then the reduced evolution of ps is not 
exactly unitary and not all the energy change can be 
identified with work. To prove this, one can compute 
the von Neumann entropy of both sides of Eq. (17) and 
notice that the mutual information I[psa) = 0. Note, 


however, that the only way this is possible is if V (C/ave) = 
Li 0 V 2 , and further we would need Vi = C/ave- But 
this is not possible, because from the fact that I/(C/ave) 
is strictly energy-preserving we can prove that Vi and 
V 2 must both be as well. Hence, the right-hand side 
of Eq. (17) cannot be a unitary if pA contains only a 
bounded reference frame. In fact, Eq. (17) can only hold 
as a limit case of using a larger and larger coherence 
resource. In summary. Assumptions 1-3 together with 
identification of work with energy change during unitary 
processes, imply that without an unbounded reference 
the work extraction protocol from Refs. [35, 36] cannot 
extract an amount of work equal to AF(ps) from a state 
with coherence. 

The strength of this point is that we do not even need 
to require the repeatability of the protocol using the same 
ancillary system. In fact, the argument is rather general. 
In order to extract all the free energy from a state ps one 
needs to transform it into a thermal state. This cannot 
be achieved by only changing the energy spectrum of FIs, 
but also requires the rotation of the energy eigenbasis, so 
that the system is incoherent at the end of the transfor¬ 
mation. This can be performed perfectly only with the 
aid of an unbounded reference frame, because it involves 
unitaries that do not strictly conserve energy.^ 

B. Extracting work arbitrarily close to the free 
energy difference 

A key fact about the Carnot efficiency is that, despite 
being achieved only by ideal heat engines that do not ac¬ 
tually exist in Nature, we can get arbitrarily close to it 
through a sequence of real engines. In a similar spirit, we 
now construct a sequence of bounded thermal machines 
getting arbitrarily close to the coherence to work conver¬ 
sion limit set by Eq. (15). The main result of this Section 
can be summarised in a non-technical way as follows: 

Theorem 1. There exists a sequence of bounded ther¬ 
mal machines approaching the ideal coherence to work 
conversion of Eq. (15) with arbitrarily high probability of 
success and with an arbitrarily small change in the quality 
parameters. The limit case is reversible. 

As an immediate consequence of the fact that the limit 
case is reversible we have: 

Corollary 1. Eq. (15) provides the ultimate limit of co¬ 
herence to average work conversion. 

In the remaining part of this section we give more de¬ 
tails about the result above, first of all specifying the 
technical claim and then the main steps of the proof (the 
details of calculations can be found in Appendix E). We 


^ A useful point of view is also given by the theory of quantum 
reference frames and recovery maps [41, 49]. 
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consider a sequence of reference states pR that approach a 
classical reference. Consider an arbitrary reference state 
Pr. We will describe it by two parameters ((A),M) ac¬ 
cording to Definition 1. 

We will now show how to perform the protocol de¬ 
scribed in the previous section to extract from any pure 
state |'0) an amount of work per copy arbitrarily close to 
the free energy difference AF{\ip)), while succeeding with 
arbitrarily high probability and changing the quality of 
the reference only by a negligible amount. For simplic¬ 
ity, define f{x) = —x — kTh 2 {x), where h 2 {-) denotes the 
binary entropy. Theorem 1 can be now made technically 
precise as follows: 

Theorem 1’. Let pR be an arbitrary reference state de¬ 
scribed by ((A),M). In the presence of a thermal bath 
at temperature T and if M is large enough, there exists 
a protocol individually extracting from M copies of l-ip) 
[given by Eq. (6)/ an average amount of work M{W), 
with 

{W) > AF(IV^)) - /(2p(l -p)(l - (A))) - 0{M-l). 
The probability of success Psucc of the protocol is 

Psucc>[l-2p(l-p)(l-(A))]“, 

and it changes the quality parameters of the reference as 
follows: 

SM = 0, (5(A) < 2^1 -Psucc 

Before presenting the proof of this theorem, let us first 
comment on its scope. Note that the same result holds 
when a reference ((A),M) is used a number of times 
M' < M, as long as M' ^ 1 (this will be clear from the 
proof). In the case in which M' > M, we can apply the 
Theorem every M uses of the reference. The changes in 
the quality parameters will eventually sum up, but the 
Theorem gives a bound on them. Also, it will be clear 
from the proof that the failure of the protocol implies a 
destruction of the coherence properties of the reference. 

We now prove the Theorem by constructing an explicit 
variation of protocol introduced in Sec. Ill A and showing 
that it performs as stated in the theorem. Recall that by 
an energy conserving unitary we can rotate l"^) around 
the 2 axis of the Bloch sphere. Hence, without loss of 
generality, we can set (^ = 0 in Eq. (6) . We then perform 
steps (1) and (2) of the protocol described in Sec. Ill A M 
times, i.e., individually processing each of M copies of {ip) 
using a reference pR described by ((A),M). The choice 
of pr ensures that during this process the reference state 
will have no population in the ground state, and so (A) 
will stay constant. Then, the hnal state of the reference 
is described with probability pf^ by pR,i = prA\^ , 
where pi = (1 — 2p{l — p)(l — (A))). Notice that having 
access to a reference described by parameters (A) and 
M such that M(1 — (A)) —0, the probability pf^ can 
be made arbitrarily close to 1. This happens because 
by taking (A) close enough to 1, we can get arbitrarily 


close to unitary evolution of a system state \tp) to a pure, 
incoherent state |1). As the final state of the system is 
almost pure, the final joint state of the system and the 
reference factorizes and an arbitrarily small amount of 
entropy is generated (the back-reaction on the reference 
is given by the Kraus operator Ai alone). 

Next, we repump the reference M-\-sa^ times, where 
M = M{l—p), um = \/Mp{l — p) and s > 0. This guar¬ 
antees that the reference has arbitrarily small population 
in states {|0)... |M)}, so that by performing a measure¬ 
ment we can project the reference to a state p'f^ with 
support on the subspace spanned by {|*)}i>M with arbi¬ 
trarily high probability (s fixes the confidence level, see 
Appendix E for details). More precisely, after repeating 
steps (1) and (2) of the protocol M times and repumping 
as explained above, the reference is described by a state 

with probability 

Psucc > [1 - 2p(l - P)(l - (A))]^E4M1/6), (18) 

where Es(x) = erf(sa:/-\/2) and erf denotes the error 
function. The final state is given by p^, described by 
((A)",M"), where M" = M and (5(A) := (A)" - (A) is 
bounded as follows 

<5(A) < 2^1-Psucc < 2^1-pfEs(Mi/6). (19) 

Notice that by taking s large enough (but hnite) we can 
make the factor Es(M^/®) in the previous two equations 
arbitrarily close to 1, Es(M^/®) « 1. In the appropri¬ 
ately chosen limit (A) —1 and M —>• oo the quality pa¬ 
rameters of the reference state are then unchanged with 
probability 1. Let us also note that the cost We of the 
measurement described above is bounded by fcT/i 2 (psucc) 
(see Appendix E). 

We have just shown that following the procedure above 
we can guarantee repeatability with arbitrary confidence 
level. Hence, we now proceed to proving that it also 
allows for extracting an average amount of work per 
system arbitrarily close to the free energy difference 
AF{\i[)). To see this, note that after repeating the pro¬ 
tocol on M copies of \i[i) we are left with M copies of 
a state T’(pg) from Eq. (9) with q given by Eq. (11), 
q = I — 2p(l — p)(l — (A)). This state is diagonal in the 
energy eigenbasis and the average work (IF) extracted 
from it is given by AF{'D(pg)): 

(IF) = l + kT\ogZ-2p{l-p){l - (A)) 

-fcrh2(2p(l-p)(l-(A))), (20) 

where as before h 2 {‘) denotes the binary entropy. By 
choosing M large enough, we can ensure that the ex¬ 
tracted work is arbitrarily peaked around the average 
given by the above equation (M can be bounded using 
the results of [9]). This ensures that when we need to 
repump the reference, we actually have enough work to 
invest to do it. The repumping costs M + S(T^^ units of 
extracted work and the cost We of the measurement is 
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independent from M. Hence, the net gain per processed 
copy of \tp) is given by 

{W) - = {W)-0 , (21) 

where 

{W) = l^F{m-2p{l-p){l-{K)) 

-kTh2{2p{l-p){l-{A))). (22) 

Therefore, the deficit per copy scales as and by 

choosing M large enough it can be made arbitrarily 
small. Moreover, the previous equation gives us the re¬ 
lation between the quality of the reference and the av¬ 
erage extracted work, showing that {W) — /S.F{\ijj)) as 
(A) ^ 1,M ^ oo, M(1 - (A)) ^ 0. 

We conclude that it is possible, with arbitrarily large 
success probability, to extract an amount of work arbi¬ 
trarily close to the free energy change from a pure state 
with coherence in energy eigenbasis, while processing it 
individually and properly taking account of all the re¬ 
sources used, i.e., ensuring arbitrarily exact repeatabil¬ 
ity. 


V. EXTRACTING WORK WITH PERFECT 
REPEATABILITY AND BOUNDED THERMAL 
MACHINES 

In the previous section we have shown how to extract 
as work all free energy of a pure quantum state with 
coherence. However we allowed for 

1. The limit case of an unbounded thermal machine, 

(A) ^ 1. 

2. An asymptotic protocol individually processing a 
large number M of copies of the system. 

These assumptions may be too strong if the reference 
itself is a microscopic system involved in the thermody¬ 
namic processing and exclude the applicability to single¬ 
shot scenarios. What if we only want to process a small 
number of systems? This requires us to go beyond the 
results of the previous section. 

Moreover, even if we only want to release the first of 
the two assumptions, i.e., put a bound on the coherence 
properties of the reference, we are still left with open 
questions. In this case the general result stated by The¬ 
orem 1 is applicable, however the work extraction proto¬ 
cols presented always entail a failure probability 1 — Psucc 
that can lead to a complete destruction of the coherent 
properties of the reference. Even if this probability is 
relatively small, we may not be willing to take this risk. 
Also, the reference inevitably deteriorates, even if by a 
small amount bounded by Eq. (19). A crucial question is 
then: are there work extraction protocols with (A) < I 
such that (5(A) = 0 and Psucc = 1? In other words, can 
we extract work from coherence using a protocol that 


never fails and gives back the thermal machine with ex- 
aetly the same quality parameters, even if the reference 
is bounded? 

In this section we construct such protocols for both av¬ 
erage and single-shot work extraction. These ensure per¬ 
fect repeatability, but the price we pay is that the average 
amount of extracted work is strictly smaller than the free 
energy difference and it is only possible for (A) above 
a certain threshold value A^rit- In the case of single¬ 
shot work extraction we show similarly that there exists 
a threshold over which the reference allows us to out¬ 
perform the single-shot protocol with no coherence. For 
clarity of the discussion, we focus on the paradigmatic 
case of the class of states I 7 ) introduced in Eq. (7). 


A. Average work extraction 

In absence of an external source of coherence no work 
can be extracted from the state I 7 ) on average [ 10 , 28]. 
However, if we allow for a repeatable use of the refer¬ 
ence, positive work yield can be obtained. In order to 
achieve this, during step ( 2 ) of the protocol we perform 
average work extraction from the state 'D{p'g) specified 
by Eq. (9). As 'F>{p'g) is diagonal in the energy eigenbasis, 
the results of [9, 10] apply. Therefore, the average work 
yield is given by the free energy difference AF{'D{pg)). 
To ensure perfect repeatability we repump the reference 
at each run, so that if (0|/9 _r|0 ) = 0, then (Ojp'^jO) = 0, 
and the reference quality parameters do not change. The 
repumping requires a unit of work, so that the work ex¬ 
tracted on average during one run of the protocol is 

{W) = q + kT{logZs-h 2 {q))-l. (23) 

The connection between the properties of the reference 
and the work yield is given by Eq. (23) together with 
Eq. (11) (where i?oo = 0 and p = r). 

In Fig. 3 we show how much work {W) can be unlocked 
through our protocol as a function of the quality of the 
reference (A) and the thermal occupation r of the excited 
state. The graph shows that the quality of the reference 
needs to be above a certain threshold in order to get 
positive average work yield. As expected, the advantage 
is the most significant for high r, because the states [ 7 ) 
and 7 s differ most in this case or, in other words, the 
amount of coherence to be unlocked is higher. 

As already mentioned in Sec. H, in the asymptotic 
regime of individually processing large number of copies 
of I 7 ), the fluctuations in the work yield, Eq. (23), be¬ 
come negligible. Notice, however, that even if {W) > 0 
we may not be able to perform step (3) every time, as 
the fluctuations around the average mean that we will 
not always have enough work to invest in the repump¬ 
ing. To resolve this problem we can follow a strategy 
analogous to the case of unbounded reference. That is, 
we repump after having extracted work M times, where 
M is sufficiently large to neglect the fluctuations around 
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FIG. 3. Coherence boost to average work extraction 
from 1 7 ). Work can be unlocked from the coherence of the 
system using a thermal machine that never deteriorates. The 
quality of the machine, measured by (A), must be bigger than 
some threshold value (boundary of the grey region) to ensure 
lyV) > 0 (r denotes the thermal occupation of the excited 
state, r = (1 + Over the threshold, the higher the 

quality the greater is the average work yield from quantum 
coherence. 

(IF).® The protocol will be repeatable up to an arbi¬ 
trarily small probability of failure, if the support of the 
reference initially starts high enough in the energy lad¬ 
der. It is important to stress, however, that the “failure” 
in this case does not entail a destruction of the coher¬ 
ence properties of the reference, as in Sec. IV B. It only 
requires the investment of extra work in order to ensure 
perfect repeatability. 

B. Single-shot work extraction 

Finally, we proceed to a fully quantum and single-shot 
protocol for an individual quantum state. This version of 
the protocol does not assume possessing an unbounded 
reference nor it requires asymptotic number of runs. In 
absence of an external source of coherence we can per¬ 
form e-deterministic work extraction from jq) - as jq) is 
indistinguishable from 75 , the results of [9] apply. This 
means that we can extract kT log Zs work with failure 
probability r or 1 -|- kT log Zs with failure probability 
1 — r. We now show that exploiting the reference in a 
perfectly repeatable way the failure probability for ex¬ 
tracting kT log Zs can be decreased - and the higher the 
quality of the reference, the stronger the improvement. 

During step (2) of the protocol we perform 
e-deterministic work extraction from the state ^{p'g) 
specified by Eq. (9), in accordance with Assumption 3. 


One can think of alternative protocols as well, in which at every 
repetition we toss a coin to decide if we repump the reference 
or not. We do not delve into this, but we expect to find similar 
results. 


FIG. 4. Coherence boost to single-shot work extrac¬ 
tion from |q). A thermal machine, used in a repeatable way, 
can exploit the quantum coherence of the system to decrease 
the failure probability in single-shot work extraction. Sim¬ 
ilarly to the case of average work extraction, the quality of 
the machine (A) must be over some threshold value to lead 
to any improvement (r denotes the thermal occupation of an 
excited state). As (A) increases, the failure probability de¬ 
creases from r to r — 5e, down to zero, i.e., to the point when 
the single-shot work extraction from the pure quantum state 
|q) becomes fully deterministic. 


With probability q we extract 1 -I- kT log Zs work and 
with probability 1 — q our protocol fails. As we need one 
unit of work to repump the reference [see Eq. (10)], the 
net gain is kT log Zs- When the protocol fails (with prob¬ 
ability q), the reference is returned in the state p'p^ and 
one has to invest one unit of work to ensure repeatability. 

In Fig. 4 we present the decrease 6e in the failure prob¬ 
ability e achieved by our protocol as compared to work 
extraction from 75 . We see that if the quality of the ref¬ 
erence is high enough, the coherence content of jq) can be 
exploited to provide an advantage in the work extraction. 
In the limit of a very high quality (unbounded) reference, 
(A) —1, the failure probability can be sent to zero, i.e., 
the work extraction from jq) becomes deterministic. 


VI. CONCLUSIONS 

In this work we have addressed the following question: 
how much work can be extracted from a state that is a 
superposition of energy eigenstates? We argued that this 
question, within the currently developed theory of ther¬ 
modynamics of individual quantum systems, is a subtle 
issue. We showed that the optimal coherence to work 
conversion can be obtained only in the limit of accessing 
a reference system with unbounded coherence resources. 
Although no real reference is unbounded (in the same 
way in which no heat engines is ideal), we can get ar¬ 
bitrarily close to the limit by means of a sequence of 
bounded thermal machines. 

The access to arbitrarily large resources should be 
questioned in the regime under study. Generally speak- 
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ing, recovering traditional thermodynamical results re¬ 
quires extra assumptions, which all entail some notion 
of “classicality”, so to effectively make coherence negli¬ 
gible: neglecting the energy fluctuations due to super¬ 
position; assuming the existence of a source of coherence 
that experiences essentially no back-reaction; collectively 
operating on infinitely many copies of the system. 

When the “classicality” assumptions are dropped one 
after the other, the results are quantitatively different 
from the thermodynamics of incoherent systems. Never¬ 
theless, we find that the coherence between energy lev¬ 
els can still enhance the performance of work extraction 
protocols. There exist perfectly repeatable processes ex¬ 
tracting on average a larger amount of work that could 
be extracted in the absence of coherence; and single-shot 
protocols in which coherence improves the success prob¬ 
ability of work extraction. Although these protocols are 
better than the correspondent incoherent ones, they do 
not achieve the performance reached in the classical limit. 

We also point out that while dealing with microscopic 
systems, the accounting of all the resources involved in 
thermodynamic processes becomes a crucial and non¬ 
trivial task. In this regard, we underline the importance 
of accounting for the resources which make up a ther¬ 
mal machine, and the concept of repeatability, that es¬ 
sentially captures the idea of using these extra resources 
without degrading them. In particular, the considera¬ 
tions here suggest that a full theory of thermodynamics 
in the quantum regime will require a better understand¬ 
ing of the accounting of coherence resources, including 
those found in the thermal machine. Some laws which 
place restrictions on coherence have been introduced in 
[ 8 , 28, 30, 32, 33], but we are still far from having a full 
understanding. We hope to have convinced the reader 
that the question of the role of quantum coherence in 
thermodynamic considerations does not admit an easy 
and immediate answer, and that it is only by appropri¬ 
ately incorporating it into the theoretical framework that 
we can explore truly quantum mechanical effects. 
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Appendix A - Collective processing regime 

In the collective processing regime work can be effec¬ 
tively unlocked from coherence. This is achieved by pro¬ 
cessing many copies of a system state ps collectively and 
extracting work from relational degrees of freedom that 
live in decoherence-free subspaces [10, 26, 28, 41]. The 


intuitive explanation is that one copy of a state ps with 
coherence can act as a reference for the other one, and we 
have ^ 'D{ps)^'^■ In the case of finite number of 

copies p®^ a non-zero amount of work is unlocked from 
the coherences, and in the limit of processing collectively 
infinitely many independent and identically distributed 
(i.i.d.) copies, the amount of work per copy that can be 
extracted deterministically equals F{ps) — F{js)- 

Instead of collectively processing many copies of a sys¬ 
tem, one may also consider a black-box device B that 
takes in individual quantum systems ps and at each 
round returns a thermalized state ys and an average 
amount of work equal to F{ps) — F{'ys)- Erom outside 
the box it seems we are dealing with a work extraction 
protocol that individually processes each state. However, 
the devious way in which the box achieves this is the fol¬ 
lowing: 

1. The box B contains a large quantum memory con¬ 
sisting of A ^ 1 copies of incoherent quantum 
states erf for which F{aB) = F{ps). 

2. Every time the box takes in a single copy ps it 
swaps this state into memory and instead performs 
work extraction on one copy oi a b- Hence it out¬ 
puts on average F{ps) — F{js) and a thermalised 
state 7 s. 

3. After N uses its memory is Hlled with the coher¬ 
ent states pg^ and so it does large-A collective 
processing and restores to erfwith costs growing 
only sublinearly with N [26]. 

Although from the outside of the box this is identical 
to the individual processing regime, the collective, rela¬ 
tional processing of coherence is “hidden” in the quantum 
memory. 

Appendix B - Average energy conservation does not 
explicitly model energy fluctuations 

Consider a system and an ancilla described by Hamil¬ 
tonians FIs and Ha, and prepared in states ps and pA, 
respectively. Assume also that the initial state of the 
system ps has coherence between energy eigenspaces. 
Now consider a joint energy-conserving unitary U, i.e., 
[U, Hs + Ha] = 0, inducing the following evolution: 

U {ps <8) paW"^ = p's® p'a^ 

so that the final state of the system pf has no coherence 
in the energy eigenbasis and pf = VpsV"^ for some uni¬ 
tary V that conserves average energy. The uncertainty 
of an energy measurement on p'a can be decomposed as 
[50]: 

H{p'a) = S{p'a) + A{p'a), (24) 

where A{a) = S{a\\'D{a)) is the relative entropy between 
a state and its decohered version and H is the Shannon 
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entropy of the probability distribution of an energy mea¬ 
surement. Because U commutes with the total Hamilto¬ 
nian we have 

A{ps® pa) = A{p's® p'a). 

As the final state of the system p'g has no coherence we 
have A[p'g ® p'^) = A{p'^). Using A{ps 0 pa) > A{pa), 
one gets that A{p'^) > Alp a)- From the invariance of the 
von Neumann entropy under unitary transformations, 
^(p'a) = ^(pa)- So we conclude from Eq. (24) 

H{p'a) > S{pa) + A{pa) = H(pa). 

Appendix C - Details of the repumping stage 

Although one could question the repumping stage de¬ 
scribed by Eq. (10), given that A is not a unitary, we 
note that this is actually not a problem. This is because 
such operation can be realized through a joint energy- 
conserving unitary between a weight system in a state 
|1) and the reference in a state The unitary is given 
by V{U) in Eq. (8), where we take U = X, the Pauli X 
operator. Then 

V{X) = a_^A + a+^A^ + |0)(0| ® |0)(0|, 

where cr+ = |1)(0| and cr_ = |0)(1|. As the reference has 
no population in the ground state, the final state of the 
weight system is | 0 ) and the final state of the reference 
is given by Eq. (10). 

Appendix D - Free energy change of the reference 


where |i) are eigenstates of the reference Hamiltonian. 
We impose the requirements Mp ^ 1 and M{l—p) ^ 1, 
where p is fixed by Eq. ( 6 ). 

We now compute the probability of the occurrence 
of the Kraus Ai on a generic reference state a. From 
Eqs. (12a)-(12b) and the fact that AA^ = I we obtain 

Pi(cr) := Tr (^Ai(tA\^ = 1 - 2 p(l - p){l - Tr (Act)), 

where recall that A = (A -|- A^)/2. Define the state of 
the reference after performing work extraction on n > 1 
qubits through the following recurrence formula 

p'S'--=Aj;-»Al + A,p'i‘-'A[. ( 25 ) 

where = p^. Because pR has initially no support in 
the first M energy levels, we can extract work from M 
qubits before there is any overlap with the ground state. 
In other words, (0|p^^|0) =0, Vn € {1,..,M}. The 
previous formula, together with Eq. (14), implies that A 
is conserved throughout the protocol. Hence we deduce 
that 

Tr [a^prAX^) = pflpR) = (1 - 2 p(l -p)(l - (A)))“. 

For notational convenience, we will now drop the explicit 
dependence of pi on pR (initial state of the reference). 
Using Eq. (25) we have 

P^R^ = PiPRA + (1 - Ipr)^ (26) 

where contains all strings of Aq’s and Ai’s different 
from the string consisting only of Ai’s and 


Denote by AFA „ the change in the free energy of the 
reference at the n-th repetition of the protocol. The total 
free energy change of the reference after M repetitions of 
the protocol satisfies 

M 

^ AF^,„ < F(p^) - F(7fl) VM, 

n—1 

where y/j is the thermal state of the reference. Hence, 
the average change in the free energy of the reference as 
M —>■ oo is 


M 


AFr := lim — 

M—>-oo Af 


afr^u= 0 . 


Appendix E - Details of approaching free energy 
limit 

We provide here the details of the repumping protocol. 
We start from a generic reference state pR such that 


Pr,i = AX^ prA\^ /p^^ . 

We can now compute A^^ ■. 

We see that is binomially distributed in the num¬ 
ber of lowering operations Ab The average number of 
lowerings is M = M{1 — p) and the standard deviation 
is (Tm = a/ Mp{l — p). We can perform a number of re- 
pumpings as in Eq. (10) as detailed in Appendix C. Let 
us denote this operation by V. We have chosen M suf¬ 
ficiently large so that the conhdence levels associated to 
gm are approximately gaussian. Hence, we can repump 
the reference M + times, which guarantees that 

the reference has arbitrarily small population in states 
|0)... \M) with a confidence level controlled by s > 0 
and increasing with M. More precisely, if Pm is the 
projector on the subspace spanned by {|0),..., \M)} and 

= I ~ Pm, 


supp(pr) nspan{|0),..., |M)} = 0, 


Tr {P^VipR,i)) > eAlsM^/yV2) := E,(M1/6), (27) 
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where erf(x) = ^ Jq exp(—t^)dt denotes the error func¬ 
tion. Now, using Eq. (26) and Eq. (27) 

:= Tr > pf E,(Mi/®). (28) 


Define 


Pr ■— 


PMnp‘'r)p^ 

Pi,Pip^R^)Ph 


This implies that in performing the two-outcome mea¬ 
surement {Pm, Pm} we would find the outcome P^ with 
probability given by Eq. (28). 

Performing such a measurement guarantees that the 
final state of the reference will have no support on a sub¬ 
space spanned by {|0 )...|M)}, similarly to the initial 
state. However, performing a selective measurement has 
a thermodynamic cost that we have to take into account. 
More precisely, such a measurement can be performed 
using an ancillary memory qubit system A described by 
trivial Hamiltonian Hj^ = 0. Then, taking the initial 
state of H to be a pure state |0), we can perform oper¬ 
ation on the joint reference-ancillary state described by 
the Kraus operators Mi = P^ (g) I and M 2 = Pm 0 cr^- 
This operation is energy conserving, as the Kraus op¬ 
erators commute with the total Hamiltonian Hfj -|- Ha- 
Hence, it is free of thermodynamic cost. Now the projec¬ 
tive measurement on states |0) and |1) can be performed 
on the ancillary memory system. Observing the result 
0 will project the reference on a subspace P^, whereas 
observing the result 1 will project the reference on Pm- 
The thermodynamic cost associated with this projective 
measurement is the cost of erasing the memory system 
afterwards. This is given by 


Now, using the gentle measurement lemma [51, 52], 
Eq. (28) also implies 


Pr-P{Pr) <2Vl-ftucc. 


(29) 


From Eq. (29), and the following characterization of the 
trace norm (see [53]) 

Up — all = max Tr {Aip — a)) 

0<A<I 


we find that 

Tr (Ap") > Tr (AiP(p^"^)) - 2v'l-Psucc 
= Tr (Apn) - 2a/1 -psucc, 

where the last equality comes from the fact that (A) is 
conserved in the protocol, up to the measurement. The 
last equation can be rewritten as 

Tr (A(pfl - p)^)) < 2 a/ 1 - Psucc- 

Exchanging the roles of pR and p'^ and introducing 


We = fcT/i2(Psucc), 


5(A) = Tr (Ap" ) - Tr (Ap^) , 


which can be made arbitrarily small as Psucc 1- Notice 
that we only needed to use a classical memory to record 
the measurement outcome, which is not in contrast with 
assumption 3 of Section HI. Also note that this cost has 
to be paid only after extracting work from M copies, 
hence the cost per copy scales as M~^ - 


we conclude 

|5(A)| < 2 -/I - Psucc- 

Using Eq. (28), this bounds the maximum allowed change 
of the quality parameter of the reference. 
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